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ABSTRACT 

We systematically reanalyze two previous observations of the black hole (BH) GX 339-4 in the very high 
and intermediate state taken with XMM-Newton and Suzaku. We utilize up-to-date data reduction procedures 
and implement the recently developed, self-consistent model for X-ray reflection and relativistic ray tracing, 
RELXILL. In the very high and intermediate state, the rate of accretion is high and thus the disk remains close 
to the innermost stable circular orbit (ISCO). We require a common spin parameter and inclination when fitting 
the two observations since these parameters should remain constant across all states. This allows for the most 
accurate determination of the spin parameter of this galactic black hole binary from fitting the Ee Ka emission 
line and provides a chance to test previous estimates. We And GX 339-4 to be consistent with a near maximally 
spinning black hole with a spin parameter a, > 0.97 with an inclination of 36 ± 4 degrees. This spin value 
is consistent with previous high estimates for this object. Eurther, if the inner disk is aligned with the binary 
inclination, this modest inclination returns a high black hole mass, but they need not be aligned. Additionally, 
we explore how the spin is correlated with the power of the jet emitted but And no correlation between the two. 


1. INTRODUCTION 


GX 339-4 is an important source for accretion physics. 
It is a recurrent low mass X-ray binary system (LMXB) 
that has a prominent disk reflection spectrum. This galac¬ 
tic black hole binary (GBHB) has undergone four outbursts 
over the last twelve years, th us making it the one of the 
most active transient systems (iPlant et al.ll2014ll . Through¬ 
out its outbursts, GX 339-4 has been observed in the very 
high/soft state (VHS), intermediate state, and the low/hard 
state (LHS) followed by a very fai nt state, but has no t yet 
been observed in the quie scent state (iHvnes et al.ll2003h . See 
McClintock & Remill^dl for a review on states. The 

mass of the BH and inclination are determined through in¬ 
teractions with the companion sta r. Through the anal ysis of 
emission lines during an outburst, [Hvnes et aH (120031) deter¬ 
mined a mass function of 5.8 ± 0.5 Mq, where the mass func¬ 
tion can be expressed in the form of 


f{M) 


Mx sin 


( 1 ) 


Mx being the mass of the compact object, q the mass ratio 
of the companion to the compact object, and i the inclina¬ 
tion. The inclination of the BH is still under some specula¬ 
tion since it may deviate from the inclination of the system 
itself, which is discussed in further detail in §3. The distance 
was determined using the optical/IR absorption along the line 
of sight to GX 339 -4. It was found to lie between 6-15 kpc 
(iHvnes et al.ll2004]) though this value is generally quoted as a 
distance of 8 knc (fzdziarski et al.ll2004l) . 

GX 339-4 is an ideal source for applying new models that 
use the method of fitting Ee Ka emission to estimate spin be¬ 
cause of its prominent reflection features. The advantage of 
using the iron emission line is that it eliminates the need to 
know the mass or distance to the object. The determination 
of mass and distance is complicated by large uncertainties for 
most objects. 
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Models that allowed for the spin of a BH to be a variable 
parameter that could be constrained directly through the fit¬ 
ting of spectra emitted from an object were not available un¬ 
til 2004. However, line models did give important early in - 
dications of very high spin parameters (iMiller et al.l 120021) . 
Since then, there have been numerous additive and multi¬ 
plicative models developed to describe the different features 
that appear in the high e nergy X-ray band. Thes e include 
model s such as RELLINE dP auser et ffll2010l 1201 3h and XIL- 
LVER (iGarcfa & Kallmanll2010ll2013h among others. The lat¬ 
est of these models to be deyeloped is the combination of 
the relatiyistic ray tracing kernel, RELLINE, and X-ray reflec- 
tion code, XILLVER , into one self consistent model RELXILL 
(IGarcfa et al.ll2014l) . 

GX 339-4 was estimated to harbor a near maximally 
spinnin g BH during an outburst in 2002-2003 (IMiller et al.l 
l2004al) . Obseryations were taken by Chandra, XMM-Newton, 
and the Rossi X-Ray Timing Explorer (RXTE). This ob- 
seryation was among the first to use burst mode aboard 
XMM-Newton, and som e aspects of the re duction and sys- 
tematics were unknown. iKirsch et al.l (l2006l) later extensiyely 
went through analyzing burst mode data taken of the Crab 
Nebula in order to al low for a more conc lusiye extraction of 
spectra. In any case, IMiller et'd] (l2004al) ruled out the pos¬ 
sibility of this object being a non-spinning BH at more than 
the 8 a confidence leyel from the iron (Ee K^) emission line 
profile. They were able to estimate that the accretion disk ex- 
tended in t o (2 — 2>)Rg, which translates to a* > 0.8 — 0.9. 
iReis et aO (120081) r eanalyzed the obseryation s of GX 339-4 in 
the VHS and LHS (IMiller et al.ll2004^ |2006^ which were ob- 
seryed by XMM-Newton. They used updated reflection mod¬ 
els (reflionx and refhiden) that were able to take into 
account blackbody emission and Compton scattering in the 
accretion disk that were not ayailable prior. The results of 
their extensiye modeling was a spin of 0.935 ± 0.001 which 
supported the preyious findings. 

These results were strengthened further by obseryations 
taken by Suzaku af t er an o utburst that occurred during 2006- 
2007. IMiller et'SH (120081) fit the spectra obtained during the 
outburst with two other spectra taken by XMM while GX 339- 
4 was in a higher and lower state. The joint fit across states 
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returned a BH with a spin parameter of a* = 0.93 ± 0.01. 
This was again among the hrst observations in a burst mode, 
and the joint fitting procedures served as an important check . 
Later analysis of the Suzaku data by lYamada et alJ (l2009l) 
found that the spectra may b e piled up; howeve r, a compre¬ 
hensive study of pile-up by iMiller et^ (120101) shows that 
pile-up only acts to narrow reflection features, and cause spin 
to be under-estimated. 

Spin m ay be hard to determ ine in the low/hard state as can 
be seen in iDauser et al.l ( 120141) and iFabian et alJ (120141) . there¬ 
fore making it necessary to focus on the very high and in¬ 
termediate states using new models to deliver the best possi¬ 
ble spin measurement for GX 339-4. This focus on the VHS 
and intermediate state has to do with the location of the inner 
disk in relation to the ISCO. This is discussed further in §3. 
We reanalyze the obseryations from IMiller et al.ll2004al and 
IMiller et^l2008l using the latest data reduction procedures 
and model. This paper is formatted as follows: §2 describes 
the obseryations and data reduction process. §3 is the spectral 
analysis and results. §4 presents a discussion of the findings 
of this work. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. XMM-Newton 

We reanalyze the 75.6 ks obseryation of GX 339-4 in 
the yery high state during reyolution 0514 that began on 
2002 September 29 at 09:06:42 UT. The EPIC-pn camera 
(IStrtider et alJl200n) was operated in burst mode with a thin 
optical Alter. The data were reduced using XMM-Newton 
Scientific Analysis System (SAS yer. 13.5.0). We also fol- 
low the gui d elines for data reduction in burst mode from 
iKirsch et alj (|2006l) . The data were reduced two different 
ways using the task EPFAST and again without the task since 
the use of EPFAST m ay alter the energy scale of the Fe K band 
d Walton et alJl201^ . RAW space was set to RAWY< 140 
and RAWX was 10 columns across (30-40). Aboye RAWY 
of 140 the determination of normalization becomes increas¬ 
ingly more affected based upon examination of figure two in 
Kirsch et. al. (2006). Although photon indices don’t become 
affected until aboye RAWY> 160, choosing to lower the limit 
ensures that F is unaffected by pile up and a better determina¬ 
tion of the normalization. In addition, the eyents were Altered 
to allow single and double eyents (PATTERN< 4) while ex¬ 
cluding bad pixels (FLAG==0). 

Since the object is so bright, there was not a region ayail- 
able to extract a background. Most obseryations of GX 339- 
4 with XMM-Newton haye operated the EPlC-pn camera in 
“timing” mode, eyen when the source is not bright, serying 
to complicate background estimates. An entire Epic-pn CCD 
is exposed in timing mode. To estimate an upper limit on 
the background, we extracted an area equiyalent to a Epic- 
pn CCD (4.4 by 13.6 arcmin) from the MOS-1 exposure of 
GX 339-4 obtained in obsID 0112900201. We chose a re¬ 
gion next to GX 339-4, and with a space density of back¬ 
ground sources that is higher than ayerage, as a kind of worst- 
case scenario. Extracting the data from this region, includ¬ 
ing point sources and any diffuse emission, giyes a flux of 
f = 1.1 X 10“^^ ergs cm“^ s“^. In fact, this limit is 
much too high, as the extraction region that we used within 
timing mode is smaller than the full CCD. Eyen so, this 
flux is four orders of magnitude below the flux measured 
(/ = 1.1576 X 10“® ergs cm“^ s“^) from GX 339-4 in the 
obseryation of interest. 


Response files were generated using SAS RMFGEN and AR- 
FGEN commands. There are at least 20 counts per energy bin 
to allow the use of minimization when performing spec¬ 
tral fits. GX 339-4 is determined to be in the yery high state 
at the ti me of this obseryatio n due to the spectral parameters 
fo und in IM iller et aH (l2004al) and the timing properties found 
in iHoman et al.l (l2005al) . 

2.2. RXTE 

RXTE simultaneously obseryed GX 339-4 on 2002 Septem¬ 
ber 29 at 09:12:11 UT for 9.6 ks. There was not a simul¬ 
taneous obseryation taken at the time of the Suzaku obser¬ 
yation, but rather, there were four other obseryations taken 
on the same day for 3.5, 3.41, 3.4, and 1.8 ks. We chose 
to use the standard products that are generated from standard 
modes after each obseryation as a guide for the photon in¬ 
dex of the XMM-Newton and Suzaku obseryations. We added 
0.6% systematic error to the spectra using the HEASOFT (yer. 
6.15) sub package FTOOL GRPPHA task. The Proportional 
Counter Array (PCA; Llahoda et al.ll200(j) spectrum is gener¬ 
ated from the standard 2 mode which coyers 129 energy chan¬ 
nels and has a 16s time resolution. The two High Energy 
X-ray Timing Experiment (FlEXTE: [Rothschild et al.l f l 99^ 
spectra (from Cluster 0 and Cluster 1) haye the same time res¬ 
olution as the PCA and coyer 129 spectral channels. Refer to 
The RXTE Cook Book for all these details. When performing 
spectral fits, the PCA was At in the 3.0-25.0 keV energy range 
and HEXTE At the 20.0-100.0 keV energy range. 

2.3. Suzaku 

We reanalyze the obseryation of GX 339-4 in the in¬ 
termediate state performed by Suzaku on 2007 February 
12 beginning at 05:33:31 UT. This obseryation is classi¬ 
fied as being in the intermediate state due to many prop¬ 
erties that are consiste nt with preyious re sults obtained in 
the intermediate states (IMiller et al.l 12004^ . We processed 
the unflltered eyent files through AEPIPELINE to produce 
clean eyent files for both the X-ray Imaging Spectrometer 
('XIS:[ifo yama et alJl200'^ and Hard X-ray Tele scope (HXD; 
iKokubun et akl 1200'i^ nd iTakahashi et akl feoOTh instruments. 
XIS was operated in 0.3s burst mode with the 1/4 win¬ 
dow option selected. HXD cameras were operated in the 
standard mode. All proper calibration products were deter¬ 
mined and applied from the HEASARC Calibration Database 
(HXD20110913, XIS20120209, and XRT20110630). 

We chose to use the two working ‘front-illuminated’ cam¬ 
eras, XISO and XIS3, in our analysis. The ‘back-illuminated’ 
instrumen t, XIS 1, wa s not u sed due to calibration uncertain¬ 
ties as per lPlant et al.1 (120151) . New attitude files were created 
for both editing modes on each camera using AEATTCOR.SL0 
to account for the wobbling that Suzaku undergoes. The at¬ 
titude files were then applied to the clean eyent files by us¬ 
ing the FTOOL XISCOORD. Eyen though the 1/4 window 
and burst option were selected for this obseryation, there is 
still a significant amount of photon pileup in the data. Pileup 
causes the photon index of a spectrum to appear harder due 
to the inability of the detector to register indiyidual eyents. 
Two or more eyents will be interpreted as a single eyent with 
a higher energy (see IMiller et'm 120101 for a full description 
of the problem and consequences). Therefore, we use the 
RXTE obseryations that were taken on the same day as a 

^ littp://space.mit.edu/ASC/software/suzalcu/aeatt.litml 
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guide for the photon index of GX 339-4 in the 2.0-10.0 keV 
energy range. We find that excluding data affected by pileup 
at the 8% level produces the proper spectral slope. We employ 
the tool PILEEST to determine the contours for this region for 
XISO and XIS3. We then exclude this region when using XIS 
XSELECT by using an annulus when extracting spectra. We 
used an outer radii of 7.5' with an inner radius of 1.1' and 1.4' 
for XISO and XIS3 respectively. There were not available re¬ 
gions in the images to extract a background. Response files 
were created for each spectrum using the task XISRESP before 
combining the two XIS spectra together with the FTOOL AD- 
DASCASPEC. Last, the combined spectrum was run through 
GRPPHA to require a minimum of 20 counts per energy bin 
like the XMM-Newton spectra. 

The HXD/PIN was operated using the XIS nominal point¬ 
ing. The proper response and non-X-ray background (NXB) 
files were obtained before executing the FTOOL HXD- 
PINXBPI on the cleaned event. This tool creates a source 
spectrum that is corrected for dead time as well as a com¬ 
bined background (cosmic and non-X-ray) spectrum. A min¬ 
imum of 20 counts per energy bin was required again. The 
GSO cleaned event file was run through the FTOOL HXDG- 
SOXBPI after the appropriate NXB, spectral binning and addi¬ 
tional ARE files were downloaded. The GSO NXB is already 
binned so the source spectrum must be binned accordingly. 
When performing spectral fits the XIS, PIN, and GSO were 
fit in the 1.0-10.0 keV, 15.0-60.0 keV, and 70.0-600.0 keV en¬ 
ergy ranges respectively. 

3. SPECTRAL ANALYSIS AND RESULTS 

We use XSPEC version 12.8.1 (lArnaudI 119961) in this 
work and errors are quoted at at least the 90 % confidence 
level. We utilize th e new RELXILL vO.le modeling package 
dGarcfa et al.l2014l) to properly describe the reflection and rel¬ 
ativistic effects within the spectra, specifically to fit the Fe Ka 
emission line and determine the spin of this object. We com¬ 
bine RELXILL with the additive model DISKBB to account 
for the thermal component of the disk and the multiplica¬ 
tive model TBABS to account for absorption along the line of 
sight. The absorption column is fixed throughout this analysis 
to be 5.26 x 10^^ cm~^ in accordance with the findings by 
[Dickey & LockmanI (fl9^ . Note that GX 339-4 was in the 
“very high” state during the XMM-Newton observation and in 
the intermediate state for the Suzaku observation. 

The inclination of GX 339-4 has not been determined con¬ 
clusively. It is known that the orbital inclination has to be less 
than 60° from the la ck of eclipses present in the optical data 
jCowlev et al.ll200^ . Luther, it can not have a inclination of 
much less than ^ 40° in ord er to have a dynami cal mass con¬ 
sistent with the findings of iHvnes et al.l (l2003h . Any lower 
of an inclination would gi ve a BH with a mass greater than 20 
Mf:^. lYamada et al.l (I2009h found that the whole 25°-45° range 
is allowed at the 90% confidence level. On the other hand, 
the orbital inclination may not even be aligne d with the incli¬ 
nation of the BH itself (see lMaccaronell^02h . We choose to 
allow the inclination to be free when fitting the spectra though 
we find it to be consistent with th e 90% confidence range de¬ 
termined by lYamada et al.l (120091) . 

We fit the two observations jointly in order to obtain bet¬ 
ter constraints on the parameters of greatest interest. The 
individual spectra provide fits of comparable quality, but do 
not give equally strong constraints on their own. Tying the 
Suzaku data with the EPIC-pn data provides a better esti¬ 
mate of the physics occurring within the source. The EPIC- 



Figure 1. Simultaneous fit of Suzaku XIS (black) and XMM-Newton (red) 
with the addition of the HXD/PIN spectra (blue) extending to the higher en¬ 
ergies. XIS data has been rebinned for plotting purposes. See Table 4 for 
parameter values. See Figure 2 for unfolded model spectrum. 



Figure 2. Unfolded model spectrum from Figure 1 for Suzaku XIS and 
HXD/PIN. The blue component corresponds to DISKBB. The red component 
illustrates RELXILL. 



Figure 3. The change in goodness-of-fit versus spin taken over 100 evenly 
spaced steps generated with XSPEC “steppar”. The spin was held constant 
at these steps while the other parameters were free to adjust. The inset shows 
a closer view of how changes for spin greater than 0.85. The blue dotted 
line in both cases is the 90% confidence level. 
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Table 1 

Joint Fitting of Suzaku and XMM-Newton with HXD/PIN 


Component 

Parameter 

Suzaku XMM-Newton 

TBABS 


0.526 


DISKBB 

TinikeV) 

0.76 ±0.01 0 

.80 ±0.01 


norm 

2230j:gJ° 

1640/^° 

RELXILL 

Qin 

6.9 

7± 1 


Qout 

U.UUUZ_Q QQQJ 

6- 

6 d 

+ 1 

CN 


^break{.^g^ 

c O + 0.4 
^•'^-2.6 

4.4 ±1.5 


a* 

> 0.97 



*(°) 

36 ±4 



R,n{iscoy 

1 



RoutiRgY 

400 



z' 

0 



V 

1.99 ±0.01 2 

.56 ±0.02 


logiO 

4.36 ± 0.04 

4 

^•'-0.3 


Xpe 

4 q +2.0 


Eaut{keV) 

+ 1 

o 

o 



freji 

2 4+0.8 

-*-•-*--0.5 


angleon' 

1 

'^•^-0.4 


norm 

0.8 ±0.2 


x/(dof) 

1.23 (1472) 


' = fixed 


Note.— Errors are quoted a t > 90 % confidence level. Nh was fixed as the 
IDickev & LockmM <199011 value. A constant was allowed to float between 
the XIS and HXD/PIN spectra. The XIS was frozen at the value of 1.0 and 
the HXD/PIN factor was fit at 0.857. A Gaussian was modeled at 1.57 keV 
with an equivalent width of 6.38 X 10“® keV and 3.08 X 10“^ 
normalization. The inner emissivity has been restricted between 3-10. The 
outer emissivity has been restricted between 0-3. The break radius was 
restricted between 3-6 Rg. qin for the Suzaku spectra is not well 
constrained within restricted boundaries. Setting angleon=l takes the 
inclination into account when modeling reflection. log{^) for XMM-Newton 
is consistent with the hard limit of 6 within the 90% confidence level. 

pn spectrum has a limited bandwidth which may not render 
the proper spectral index. In order to obtain a better esti¬ 
mate of the photon index, F, we fit a grid of RELXILL models 
(for a* = 0,0.5,0.99 and i < 45) with a blackbody disk 
component to the RXTE spectra for the simultaneous obser¬ 
vation taken in 2002. A constant was allowed to float be¬ 
tween the PCA and HEXTE. The photon index, F, was de¬ 
termined to be 2.56 )')q q 3 , 2.61 ± 0.03, and 2.621 q q 5 for the 
different values of fixed spin. To account for calibration un¬ 
certainties between the two missions, we restrict the photon 
index of the EPIC-pn spectrum to be within AF=0.1 of the 
average value of F = 2.60 when applying models. Again, 
the XMM-Newton data were reduced two different ways with 
and without the task EPFAST. We initially used the data that 
were reduced with EPFAST when modeling the data jointly 
with Suzaku. We then overlaid the spectrum that was gen¬ 
erated without EPFAST onto the best fit to see if there was 
a significant difference between the two XMM-Newton spec¬ 
tra. The differences were very small so we chose to use the 
data generated with the task EPFAST. We ignore the 1.75-2.35 
keV range in the XMM-Newton data rather than fitting mul¬ 
tiple Gaussians to accou nt for features of instrumental origin 
as per lPlant et al.l (l2015h . Additionally, we ignore below 1.3 
keV d u e to residuals that wer e also identified b y iReis et aP 
(1201 Ih . iHiemstra et al.l (1201 Ih . and iPlant et ^ (1201 5ll . For 
the Suzaku data, we ignore the 1.5-1.9 keV and 2.1-2.4 keV 
energy ranges due to the poorly calibrated known Si-K and 
Au-M instrumental residuals. 

As stated before, the model that we use to fit the spectra 


is TBABS*(DISKBB-IRELXILL). In Table 1 are the values ob¬ 
tained when performing a joint fit between the XMM-Newton 
and Suzaku observation in the l.O-lO.OkeV energy range and 
extended out into the higher energy range from 15.0-60.0 keV 
with the spectra obtained from the HXD/PIN. The spin and 
inclination have been tied between the two observations be¬ 
cause those should remain consistent between observations 
and state changes. The angleon parameter can be set to ei¬ 
ther 0 or 1. Setting the angleon parameter to 0 causes the 
model to angle-average the X-ray reflection. Using an angle- 
averaged X-r ay reflection model can bias the inclination to 
lower values dGarcfa et al.ll2014h . We make sure to set the 
angleon=l so that the inclination is properly taken into ac¬ 
count and not angle averaged. Figure 1 shows the simultane¬ 
ous modeling of the spectra from XMM-Newton and Suzaku. 
Figure 2 shows the unfolded model spectrum for the Suzaku 
XIS and HXD/PIN as reference for what the models are ac¬ 
counting for the data. 

We also assume in this model that the inner accretion disk is 
fully extended to the ISCO. Disk truncation is still a topic of 
debate, though there has been some recen t evidence presented 
for truncation in the intermediate state (iTamura et al.ll2Ql^ 
iPlant et ^12015h. it seems to be unlikely even in the low/hard 
stat^ (see iMiller et al.l 1200^ iRvkoff et al.l UOOTt iReis et aP 
l2009ll2010l) . A strong indicator that the inner disk is extended 
down to the ISCO i s the presence of high frequency quasi- 
periodic oscillations (lNowakl200C)llNespc)li et al.ll2003h in the 
intermediate state and the VHS. The freq uencies are propor- 
tional with the orbits close to the IS CO (IStrohmayenl200lL 
IMiller et aklllioTl iHoman et al.llMil l2005bl) . We Slow the 
emissiyity, q, to be broken since constant emissiyity through¬ 
out the disk is often thought to not be physically reasonable 
due to the differ ence in conditions in th e accreting material 
closer to the BH (I Wilkins & Fabianl2012^ . For this reason we 
then allowed the inner emissiyity {qin) to yary between 3-10 
Rg and the outer emissiyity index (gout) to yary between 0-3 
Rg. The break radius was restricted to lie between 3-6 Rg. 

It can be seen that the model returns an acceptable fit 
to the data (x^/dof = 1816.36/1472). The hi gh ioniza¬ 
tion pa rameter is consistent with those obtained in IReis et al.l 
(120081) when looking at obseryations taken in the VHS and 
LHS. It was found that ionization parameter for the VHS was 
I'Ogi.O > 4 and in the LHS to be log{^) « 3. It makes 
sense that the ionization parameter in the intermediate state 
would fall between the two yet be closer in yalue to the VHS 
since it was taken after the outburst had peaked. The ioniza¬ 
tion parameter yalue for the XMM-Newton spectrum is con¬ 
sistent with the hard limit of 6 within the 90% confidence 
leyel. The iron abundance was allowed to be free and returns 
a yalue that is « 4 times the solar abundance. This may be 
an oyerestimate that ca n b e attributed to the RELXILL model. 
See iGarcia et al.ll2014i and iKara et aT1l2015l for further expla¬ 
nation. The inner emissiyity for the Suzaku spectra is not 
well constrained but it is clearly not driying the spin result. 
The temperature o f the disk is consistent with those found in 
IMiller et"^ (l2008l) . The high reflection fraction suggest that 
the spectra is reflection do minated, but this is indicatiye of a 
BH w ith a high spin yalue (iDauser et al.l 12014 iParker et'dl 
12014) . The constant between the HXD/PIN and XIS is low 
with respect to the standard ~ l.l(fl but this could be from 
extracting an annular region and haying used an uncommon 
XIS mode. 

JX-ISAS-SUZAKU-MEMO-2008-06 by Maeda et a). 
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We find the inclination of GX 339-4 to be 36 ± 4 d egrees . 
Using the mass function of 5. 8 ifc 0.5 Mrr^ (iHv nes et alJ l2003h 
and a mass fraction of 0.125 (iMunoz-Darias et al.ll2008l) for 
this binary system with the inclination above suggests a large 
BH mass of ^ 37.5 Mq. This applies if and only if the incli¬ 
nation of the binary system and the BH are aligned, but this 
need not be the case. A lower mass estimate can be achieved 
if the the inclination of the BH and the binary are misaligned. 
The true inclination is hard to determine as discussed earlier 
in this section. A higher inclination would also return a lower 
BH mass. 

We obtain a spin parameter of a, > 0.97. This confirms 
the high spin estimates that have been found for GX 339-4 in 
the previous ye ars. The maxim um value for spin allowed for 
a BH is 0.998 (lThornelll974ll . When constraining the error 
on spin we find that the maximum allowed value is within the 
90% confidence level. We therefore quote the difference be¬ 
tween the maximum allowed value and value found for spin as 
the upper bound. Figure 3 shows the change in when step¬ 
ping the spin parameter from 0 to 0.998 using the “steppar” 
command in XSPEC. 

4. DISCUSSION 

We have found the spin of GX 339-4 to be a near maxi¬ 
mally spinning BH with a spin parameter of o* > 0.97. The 
high reflectio n fraction also indicates a rapidly spinning BH 
(iDauser et al.l 12014 iParker et~an |2014|). This is consistent 
previous h i gh spi n estimates made bv iMiller et alJ (l2008h and 
iReis et ahl (120081) . The analysis made in previous works and 
our analysis performed with RELXILL are different in several 
respects: 



log a. 



• REFHIDEN allows for a different atmosphere, heated by 
a strong blackbody in the mid plane. 

• RELXILL has a higher spectral resolution, improved 
atomic data, and includes more lines than prior mod¬ 
els. 

• RELXILL can be set to take the inclination angle into 
account when modeling X-ray reflection, whereas RE- 
FLIONX and REFHIDEN are angle-averaged. 

It is important to account for systematic errors. Systematics 
owing to any violation of the test particle ISCO by an actual 
fluid disk are likely imp ortant, but appear to be sm all in recent 
numerical simulations (iRevnolds & Fabi^l2008h . We have 
found that the combination of atmospheric structure, model 
resolution, and inclination angle - each potentially an impor¬ 
tant source of systematic error - only produce ^5% offsets 
in combination. In other words, they are not large sources of 
systematic error. It should be noted that inclination values that 
would keep the mass of the black hole in GX 339—4 below 
10 Mq would necessitate eclipses, which are clearly not seen. 

Given that GX339-4 is an important source, and that our 
results again point to a very high BH spin value, it is worth 
considering the role of spin on GX 339-4 and other stellar- 
mass BHs. In particular, we look into where GX 339-4 falls 
according to existing relations between BH spin and jet power 
pertaining to BH binary (BHB) systems. One such relation in 
iNaravan & McClintockI (120121) has found a near linear corre¬ 
lation for five BHBs with low mass companions and that pro¬ 
duce ballistic jets during outburst. This ap pears to confirm the 
predictions of lBlandford & ZnaiekI (1197'^ (BZ) where the rel¬ 
ativistic jets are being powered by extracting energy from the 


Figure 4. Black hole spin versus jet power. See Table 2 for values. GX 
339-4 is given in red. Errors are plotted for the spin of GX 339-4 though 
they are to small to see. The top figure is for 7 = 2. The slope of the line 
is 2.12 ± 0.75 which is close to that predicted for the BZ effect. The lower 
figure is for 7 = 5 and has a slope of 3.12 di 0.39. 

angular momentum of the BH, th ough it waslPenr ose| (|1969h 
who first theorized this concept. INaravan & McClintockI do 
not include GX 339-4 in their correlation, but provide a rough 
estimate of the spin using their correlation. They estimate GX 
339-4 to have a spin between 0.2-0.5 whe n assuming that th e 
maximum observed radio flux is 55 mJy (iGallo et al.ll2004l) . 
which differs from our findings significantly. We take the ob¬ 
served fl ux density (at 5 GHz) and spin values that are pre¬ 
sented in INaravan & McClintockI (1201 2h for four of the five 
sources: A0620-00, XTE J1550-564, GRO J1655-40, GRS 
1915H-105. These spins were determined via the continuum 
fitting method. For the maxim um observed radio flux density 
of GX 339-4 we use 55 mJy (IGallo et al.l 12004ft . Recently, 
iMorningstar & Milled (120141) conducted a spin estimate for 
4U 154747 using the iron line emission similar to what we 
have done here for GX 339-4. We choose to use this recent 
estimate when looking into the relation between spin and jet 
power. The radio flux of 4U 1543-47 is seen as a lower limit 
for the jet power and is indicated as such in Figure 4 since it 
lacks s ufficient radio d ata throughout the entirety of its out¬ 
bursts (iPark et al.ll2004l) . 

To change the observed fluxes into jet power we have to 
corre ct the fluxes for Doppler boo sting. We use an equation 
from iMirabel & Rodrfguea d 19991) to correct for this effect. 
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Table 2 

Parameters for BH spin and Jet Power 


BH Binaiy 


a* 


II 

Pieti-y = 5) 

ref 

A0620-00 

0.12 

± 

0.19 

0.12 

1.6 

1 

XTEJ1550-564 

0.34 

± 

0.24 

9.1 

220 

1 

GRO J1655-40 

0.7 

± 

0.1 

63 

650 

1 

GRS 1915+105 

0.975 

± 

0.025 

88 

1900 

1 

4U 1543-47 

0.43 

± 

0.31 

0.14 

0.53 

2 

GX 339-4 

> 

0. 

97 

0.75 

3.4 

3 


Note.— The units of Pjet are given in kpc^ GHz Jy The reference 

refers to where the spin and/o r o bserved maximum radio flux w as taken. (1) 
INaravan & McClintockll2012l (2) IMomingstar & Millei1l2014l (3) IGallo et al.l 

12004 

where the Doppler factor <5 = ( 7(1 — f3 cosO))~^, 7 = 
iV^r ^ is the Lorentz factor, /3 = v/c, and a is the 
radio spec t ral ind ex. We use updated values for a from 
iKing et aT] (l2013all . These calculations were performed for 
two different values of Lorentz factors (7 = 2 & 7 = 
5). 7 = 2 is chosen analo gous with the work done in 
INaravan & McClintockI (l2012h . but since the Lorentz factor 
is difficult to determine we choose a higher value of 7 = 5 for 
comparison. We can then use the corrected flux for each ob¬ 
ject to estimate to the power of the jet being emitted from 
the BH by using the relation given INaravan & McClintockI 

dMl). 

Pjet = D^vSv/MbH ( 3 ) 

Where D is the distance to the object, v is the observa¬ 
tion frequency (taken to be 5 GHz), Si, is the corrected radio 
flux density, and Mbh is the mass of the BH . For GX 339- 
4 we use D — 8 kpc (IZdziarski et al.ll2004ll with Mbh = 
5.8 ± 0.5Mq (iHvnes et al.l 20031) . The values for jet power 
can be seen in Table 2. Figure 4 shows the jet power in rela¬ 
tion to the spin of the BHB on a logarithmic scale. The slope 
of weighted linear regression line for 7 = 2 is 2.12 ± 0.75, 
which is roughly consistent with that predicted of the BZ ef¬ 
fect. However, the Lorentz factor is very difficult to c onstrain 
in ma ny systems and is therefore largely uncertain (iFended 
I 2 OO 3 II . When we allow the Lorentz factor to increase to 5, the 
slope of the best fit line becomes 3.12 ± 0.39. Yet, this steeper 
slope is still within 3a from the value at 7 = 2. This is not 
highly significant and can still be considered consistent with 
the BZ effect. 

It can be seen in either case that GX 339-4 lies significantly 
below the predicted jet power via spin energy extraction for a 
near maximally spinning BH. If we use the higher mass esti¬ 
mate of ~ 37.5 Mq, when assuming the the inclination of the 
binary and BH are aligned, then the source is positioned even 
further away the from the correlation between o* and Pjet ■ 
One explanation for the poor fit co uld be that the jet p ower 
is not correla t ed wi th the spin alone. iFender et alJ (1201 Ol) and 
[Russell et al.l (120131) have found that spi n and jet pow er are not 
correlated when a larger sample is used. iKing et all (2013a,b) 
have investigated the role of spin in driving jets across the BH 
mass scale. They find that spin may set the maximum power 
in the jet, but that the mass accretion rate may act as a ‘throt¬ 
tle’ in setting t he instantaneous jet pow er, as measured via 
radio flux. Both lRussell et al.l ( 1201 3l) and IKing et alJ (12013al) 
have found that a positive correlation between Pjet and a* 


is, at best, marginally significant for BHXBs (< 2.3cr confi¬ 
dence level). Our results are consistent with these findings. It 
is clear that the instantaneous jet power in GX 339-4 is not set 
by spin. 

We would like to thank the anonymous referee for the help¬ 
ful and constructive report. This research has made use of data 
and/or software provided by the High Energy Astrophysics 
Science Archive Research Center (HEASARC), which is a 
service of the Astrophysics Science Division at NASA/GSEC 
and the High Energy Astrophysics Division of the Smithso¬ 
nian Astrophysical Observatory. 
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